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Abstract 

We introduce a family of abelian sandpile models with two parameters n, m G N defined on 
finite lattices on d-dimensional torus. Sites with 2dn m ov more grains of sand are unstable 
and topple, and in each toppling m grains dissipate from the system. Because of dissipation in 
bulk, the models are well-defined on the shift-invariant lattices and the infinite-volume limit of 
systems can be taken. From the determinantal expressions, we obtain the asymptotic forms of 
the avalanche propagators and the height-(0, 0) correlations of sandpiles for large distances in the 
infinite-volume limit in any dimensions d > 2. We show that both of them decay exponentially 
with the correlation length 

^(d, a) = (v^sinh”^ a/ a{a + 2) )“^, 

if the dissipation rate a = m/{2dn) is positive. By considering a series of models with increasing 
n, we discuss the limit a J, 0 and the critical exponent defined hy Va = — liniaio log ^(d, a )/ log a 
is determined as 

I'a = 1/2 

for all d > 2. Comparison with the g 0 limit of g-state Potts model in external magnetic field 
is discussed. 

Key words. Abelian sandpile models, Dissipation, Avalanches, Height correlations, Determinantal 
expressions, Correlation length exponent. 
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1 Introduction 


Let d G {2,3,... } and L G N = {1,2,3,... }. Consider a box in the d-dimensional hypercubic 
lattice Bl = {—L, —L -|- 1,. .., C where Z denotes the collection of all integers. We impose 
periodic boundary conditions for all d directions and obtain a lattice on a torus (toroidal), which is 
denoted by A^. The number of sites in is given by |Ax,| = (2L -|- 1)'^. In the present paper we 
study a family of Markov processes on A/,, ht = {/it(z)}zgAi,, with discrete-time t G No = {0} U N. 

Assume n, m G N and let 


a = 


m 

2dn 


and 


he = 2d{l -|- a). 


Define a real symmetric matrix with size {2L 1)'^, 


he} 

AL(x,y) = <( -1, 

0 , 


if x = y, 
if |x-y| = l, 
otherwise, 


( 1 . 1 ) 


where x = (xi,... ,Xrf),y = (yi,... ,yrf) G A^ and |x - y| = “ ViY- Let l(a;) be the 

indicator function of an event a;; l(ti;) = 1, if w occurs and 1 ( 0 ;) = 0, otherwise. The configuration 
space is 

^ r 1 2 , 

= SO,he- \ . 


n n n j 

Given a configuration ht G Si^t G Nq, ht+i £ Sl is determined by the following algorithm. 


(i) Choose one site in A^ at random. Let x be the chosen site and define 

yfiAz) = ht(z)il(z = x), zgAl- 
If y^^)(x) < he, then = {r/^^^(z)}zgA^ G Sl- In this case, we set ht+i = 

(ii) If r]^y{x) = he, then r]^y ^ Sl- In this case, we consider a finite series of configurations 
{hfi)) ■ ■ ■ ’hfr)} 3r G N recursively as follows. Assume that ^ Sl with £ > 1, then 

= {z G Ai : ??^)(z) > he} / 0 and define 

yf£+i)(z) = y^^)(z) - AA(y,z), zgAa. 
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Figure 1: A toppling for the DASM with the parameters d = 2,n = 2 and m = 1. In this case 
he = 2dn + m = 9, and thus the site x with height h(x) = 10 is unstable. In a toppling, he = 9 
grains of sand drop from the site x, in which n = 2 grains land on each nearest-neighbor site, m = 1 
grain is dissipated from the system, while h(x) — he = 1 grain remains on the site x. 

If vfe+i) £ <Sl, then t = i + l and ht+i = hfr)- Remark that r = T{^,ht) and r < oo by 
Z^zizeAi ^L{y, 2 ) > 0,Vy G A^ as explained below. 

We think that 1/n is a unit of grain of sand and ht{z)n represents the height of sandpile at site 
z measured in this unit. The step (i) simulates a random deposit of a grain of sand. In the step 
(ii), for each 1 < £ <t, the sites y G A^^-^{ht) are regarded as unstable sites and the process 

{^f£)(z)}zeAi ^ - AL(y,z)} 26 Ai, 

is called a toppling of the site y such that 

AL{y,y)n = heU = 2dn + m grains of sand drop from the unstable site y 

and 

|AA(y,z)|n = n grains of sand land on each nearest-neighbor site z, |x — z| = 1. 

Since there are 2d nearest-neighbor sites of each site, m grains are annihilated in a toppling. (See 
Fig.l.) The total number of grains on Al decreases in each toppling and it guarantees r < oo. The 
conhguration space 5/, is a set of all stable conhgurations of sandpiles in which height of sandpile 
is less than the threshold value he at every site; h{z) < he, Vz G A^. From a stable conhguration ht 
to another stable conhguration ht+i, topplings occur. Such a series of toppling is 

called an avalanche. (Note that, if r = 1, toppling does not occur. Even in such a case, we call the 
transition from ht to ht+i an avalanche, which is just a random deposit of a grain of sand.) Dehne 

t(x,/i) —1 

r(x,y,h)= l(y G A^^)(h)), x,y G A^, h G Sl- (1.2) 

e=i 

This is the number of topplings at site y G A^ in an avalanche caused by a deposit of a grain of 
sand at a site x G A^ in the conhguration h G Sl- 

We have assumed that n, m G N in the above dehnition of processes. If we set n = l,m = 0, 
however, we have a = 0 and A£,|a=o gives the ‘rule matrix’ of the sandpile model introduced by Bak, 
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Tang and Wiesenfeld (BTW) [2, 3]. The BTW model have been studied on finite lattices with open 
boundary conditions in order to make r be finite. For example, the BTW model is considered on a 
box Bl. The boundary of box Bl is given by dBi = {y = ( 2 / 1 , • • • , yd) £ B^ : 1 <3i < d s.t. yi = 
—L or L}. In the BTW model defined on Bl, YIz-z&Al ^iU=o(y)2;) = 0 if y G Bl \ dB^; that 
is, the number of grains of sand is conserved in any toppling in the bulk of system. By imposing 
the open boundary condition, we have YIz-z&Kl ^iU=o(y)Z) > 0 for y G OBl and dissipation of 
grains of sand can occur in topplings at the boundary sites. In the present model, in every toppling 
at any site y G A/,, YIz zgAl ^Liy,'^)n = m grains of sand dissipate from the system and hence 
T < 00 is guaranteed in the shift-invariant system. The quantity a indicates the rate of dissipation 
in a toppling. 

The present process belongs to the class of abelian sandpile models (ASM) studied by Dhar [6]. 
We define the operators {a(x)}xeA^ following Dhar by 

ht+i = a(x)/it, X G Al, 

where ht, ht+i G Sl and the site x is the chosen site in the first step (i) of the algorithm at time t. 
That is, a(x) represents an avalanche caused by a deposit of a grain of sand at x. Then the above 
algorithm guarantees the abelian property of avalanches (see Lemma 2.1 in Section 2.1) 

[a(x),a(y)] = a(x)a(y) - a(y)a(x) = 0, Vx,yGAL. (1.3) 

We call the present Markov process the d-dimensional dissipative abelian sandpile model (DASM 
for short). The two-dimensional case was studied numerically [10] and analytically [30, 28, 18]. In 
the present paper, we will discuss the models in general dimensions d > 2 in finite and infinite 
lattices. See also [29]. As shown in [17, 26, 16] the DASM is useful to construct the infinite-volume 
limit of avalanche models. Importance of the abelian sandpile models in the extensive study of 
self-organized criticality in the statistical mechanics and related fields is discussed in [25]. 

2 Basic Properties of Dissipative Abelian Sandpile Model 

2.1 Abelian property 

First we prove the abelian property of avalanches (1.3). 

Lemma 2.1 (Dhar [6]) Assume that the avalanche operators {a(x)}xeA^ act on Sl- Then 

[a(x),a(y)] = 0, yx,y £ Al- 

Proof. Let Xl = Define three sets of maps from Xl to Xl; {t(x)}xeAi, {t(x)}xeAi and 


{d(x)}xeAi as follows, 

. For X G Ai 

and rj = 

{d(x)}xeAi 

G Xl define 


t(x)? 7 (z) = 

7?(z) - 

<1 



t(x)? 7 (z) = 

f d(z) 
1 d(z) 

<1 

1 

if ??(x) > he, 
otherwise, 


d(x)? 7 (z) = 

7?(z) -F 

-l(z = x), 
n 

z G Al- 

By definition of t. 






t(y)t(x)7?(z) 

= d(z) - 

1 

<1 

AL(y,z), z G Al 
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Similarly we have 


z G Al- 


t(x)t(y)r?(z) = r?(z) - AL(y,z) - Al(x,z), 

Therefore t(y)t(x)r/ = t(x)t(y)ry, Vr/ G that is 

[t(x),t(y)] = 0, Vx,yGAL. (2.1) 

Assume that y 7^ x. Then 

t(y),(x) = ,(x) - A(y,x) = | j* : = (• 

It implies that if t/(x) > he then t(y)T/(x) > /ic,Vy 7 ^ x, that is, any site cannot be stabilized by 
topplings which occur at other sites. Therefore, the definition of t(x) and (2.1) give 


[t(x),t(y)] = 0, 

It is obvious that 

[t(x),d(y)] = 0, 

Consider the situation that h G and 
independent of the order of the products of t(z)’s. 


Vx,y G Al. 

(2.2) 

Vx,y G Al- 

(2.3) 

7^ 0, 1 < £ < r. 

By (2.2), nz:zeA-)(h)t(z) is 


Then we can write 



r—1 

)1 

a(x)/i = 

n 

n 


e=i 

\z:zeA^^fih) J 


d(x)/i. 


X G Al, 


h G Sl- 


By (2.2) and (2.3), the lemma is proved. | 


2.2 Recurrent configurations 

Consider a subset of defined by 

TZl = {h G Sl : Vx G AL, 3 k{x) G N, s.t. (a(x))^*^^^/i = h}, 
which is called the set of recurrent configurations. 



Figure 2: The set of recurrent configurations TZl is closed under avalanches. 
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Lemma 2.2 (Dhar [6]) Ifh& TZl, then a(x)/i G TZl for any x G That is, TZl is closed under 
avalanches (see Fig.2). 

Proof. By definition, if /i G TZl, then for any y G A^, 3k{y) G N, s.t. (a(y))^(y)/i = h. If we 
operate a(x),x G Al on the both sides of this equation, then we have a(x)(a(y))^^y^h = a(x)h. By 
Lemma 2.1, LHS= (a(y))^*^^^a(x)/i. This equality implies that a(x)h G TIl- Since it is valid for 
any x G Al, the proof is completed. | 

Consider a (2L + 1)‘^-dimensional vector space Vl, in which the orthonormal basis is given by 
{e(z)}zgAi- For each conhguration r] G A^, we assign a vector 

V= Y1 (2-4) 

z:zeAx, zizSAi^ 

where 1/n denotes the unit of grain of sand. Assume that h G TZl', for each x G Al, there is 
k{x) G N such that 

(a(x))^Wh = h. (2.5) 

Consider the vector corresponding to the configuration (d(x))^(^^/i, 

T] = (h{x) + e(x) ^ h(z)e(z) G Vl- (2.6) 

^ z:z^x 

Then (2.5) claims that there exists a set {r(z) G N : z G A^} such that 

ii = V+ ^(y)^i(y>z) j e(z). (2.7) 

z-.zgAl \y:yeAi, J 

Note that (2.7) is written as 

h = ^7+ ^(y)’^(y) 

y:ysAi 

with 

v(x) = ^ AL(x,z)e(z), xgAl. (2.8) 

z:zeAi, 



Figure 3: Hypercubic lattice fl with the basis {v(x)}xgAi in Vl- Every avalanche from an unstable 
conhguration r] given by (2.6) to a recurrent conhguration h G TZl is represented by a lattice path 

77 ^ h on n. 
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We can say that, given h G TZi, all points { 77 } given by (2.6) are identified with sites of a 
hypercubic lattice Q with the basis {v(x)}xeAi in Vl. (See Fig.3.) Consider a primitive cell 
(fundamental domain) of the lattice defined by 


l^L 


< ^ c(x)v(x) : 0 < c(x) < l,x G Al > 

x:xeA/, 


cVl. 


(2.9) 


By definition, the intersection of the lattice hi and is a singleton, say p. We assume that the 
origin of this lattice is given by p and express the lattice by hl^. We consider a collection of 
all lattices with the same basis (2.8) having distinct origin in Lii, {hl^ : p G Ul}- Then there 
establishes a bijection between TZ^ = {h} and {hlP : p G 


Lemma 2.3 (Dhar [6]) The number of recurrent configuration is given by 

|7^i| = det Al. 


Proof. The above bijection implies \R.l\ = |{hlP : p G Ul}\. Since the unit of grain of sand is 1/n, 
the origins {p} of lattices {hl^} should be in (Z/n)^^, and hence 


= X (the volume of Ul). 


|{h!P :pgZYl}| = UrPifLlnt^ 

The volume of Wl given by (2.9) with (2.8) is det and the proof is completed. 



Figure 4: A primitive cell of hi on the lattice (Z/n)^^. Since the unit of grain of sand is 1/n, the 
origin p of lattice Pi should be at a site of (fLlrC)^^. 


2.3 Stationary distribution 

For /i G let be the probability law of the DASM starting from the configuration /iq = h. 

Definition 2.4 If we restrict {a(x)}xeA 2 . to IZl, inverse of the avalanche operator can be defined 
by 

a(x)'i = a(x)'^W-\ xgAl. 

Assume that h ^IZl is given. Define 
pit{X)=P^{ht = X), 

W(A ^Y)= ¥^{ht+i = Y\ht = X), t G No, X, T G Ul- 
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Consider the Master equation 


^it+i{X) = fitiX) - fit{X)W{X ^Y)+ Y1 f^t{y)W{Y ^ X), 

Y-.Y&TIl Y-.Y&Il 

where we have used the assumption that ho = h £ TZl and Lemma 2.2. By definition of the DASM, 
we can find that, for A, y £ TZl, 

W{X^Y) = Prob(x is chosen)l(a(x)A = y) 

xixsA^ 

= LO E i(=WA' = r) 

' ' x;xeAi:, 

= pTW E i(^ = a-(x)r). 

' ^ x:x&Al 

Then we have 

- fLtiX) = ^ {^i(a(x)-iA)-^i(A)}, VAG7^L. 

^ ^ x;xeAi:, 

It implies that the uniform measure on TZl, 

is a stationary distribution of the process. 

Lemma 2.5 The DASM on Kl is irreducible on TZl- 

Proof. Consider the configuration h £ Sl, such that /i(x) = he — l/n,Vx G Al- Now we take two 
arbitrary configurations X and Y from TZl- We have 

h= [J (a(x))'^^-i/’^-^WA = Yl (a(x))^=-i/’^-^Wy. (2.10) 

x:X{x)<hc — l/n x:Y{x)<hc — l/n 

Since this means that the configuration h is reachable form X and Y by avalanches, Lemma 2.2 
guarantees that h £ TZl- Since we have assumed that Y £ TZl, (a(x))^(^)y = Y with some /c(x) £ N 
for any x £ Al- Therefore, the second equality of (2.10) gives (see Definition 2.4) 

y= Yl (a(x))^W-('^=-i/'^-^W)h. (2.11) 

xiV {x)<hc. — l/n 

Combining (2.10) and (2.11) gives 

y= Yl (a[^))Y^)-ihc-l/n-Yix)) (a(y))'^--V«-^(y)x 

x-.Y{x)<hc — l/n y:X{y)<hc — l/n 

Let u = Ex:Y{x)<fec-i/n{^(x) “ (^c “ l/n - y(x))} + Ex:X(x)<Ac-i/n{^c “ l/n - A(x)}. Then we 
see 

¥'^°iht+s = Y\ht = X)> fors>a. 
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Since RHS is strictly positive for finite L, this completes the proof. | 

Then the following is concluded by the general theory of Markov chains (see, for example, 
Chapter 6.4 of [12]). 

Proposition 2.6 The stationary distribution of the DASM is uniquely given by the uniform mea¬ 
sure on TZl- 

We write the probability law of the DASM on Al in the stationary distribution as and its 
expectation as E^. 

2.4 Allowed configurations and spanning trees 

Dhar also introduced a subset of 5^ called a collection of allowed configurations Al [ 6 ]. He defined 
that for G if there is a subset F d Al such that F / 0 and 

^(y) < X] (-^i(x,y)), VyGF, (2.12) 

then h G Sl has a forbidden subconfiguration (FSC) on F. Then define 

Al = {h G Sl ■ h has no FSC}. 


Lemma 2.7 For the DASM on Al, 

TZl C Al- 

Proof. In the proof of Lemma 2.5 we have shown that h G TZl and all recurrent stares are reachable 
from this configuration h. We can prove that h G Al as follows. We assume that the contrary; 
there exists a finite nonempty set F C Al satisfying (2.12). In the DASM, however, for any 
y G F, h{y) = he — 1/n = 2d + {m — T)/n > 2d > Ylx-x&F-xj!=yi~^Li^,y)), which contradicts our 
assumption. Since both TZl and Al include h, it is enough to show that Ml is closed under the 
process of avalanche to prove the lemma, since we have already proved that TZl is so in Lemma 2 . 2 . 
Remark that addition of particles only increases h and such procedure on an allowed configurations 
cannot create any FSC. Here we assume that there exists an allowed configuration h such that by 
a single toppling at the site x it becomes to contain a FSC. Write h' = t(x)d(x)/i, that is, 

h'{y) = h{y) + -l(y = x) - Ai(x, y), Vy G Al- (2.13) 

n 

By assumption, there exists F 7 ^ 0 such that 

h'{Y)< ^ (-AL(z,y)), VyGF. (2.14) 

z:zSF:Z7^y 

Combining (2.13) and (2.14) gives 

Ky) < ^ (-AL(z,y)) + Az,(x,y), VyGF\{x}. 

z:z£F,zj^y 

Since A£,(x, y) < 0 for x / y, this inequality means that h has a FSC on F\{x} and this contradicts 
our assumption that h is allowed. Since any avalanche consists of addition of a particle and a series 
of topplings, the proof is completed. | 


9 



Definition 2.8 Given a pair {Al,Al), let U {r} with an additional vertex r 

‘root’), and he the collection 0 / |Ai(x, y)|n = n edges between x, y G A^^jX / y, 

'^yy&ts.L ~ ^dges between x G A^, and r. (See Fig.5.) Graph Gl associated to (A^,, 

is defined as 


Gl 


{dp.. 


G 


(e)', 


(the 


and 

Al) 



r 


Figure 5: A part of the graph Gl = associated to the DASM {Al,Al) is illustrated 

for the case that d = 2,n = 2 and m = 1. In this case, each pair of the nearest-neighbor vertices 
are connected by n = 2 edges and each vertex is connected to the ‘root’ r by m = 1 edge. 


Definition 2.9 We say a graph T on Gl is a spanning tree, if the number of vertices of T is 
|G^^| = |Ai| + 1, the number of connected components is one, and the number of loops is zero. 

Lemma 2.10 Let Tl = {spanning tree on Gl associated to (A/,, A^,)}. Then 

\Tl\ detAi. 

Proof. See p.l33 of [20] and Theorem 6.3 in [4]. 


Lemma 2.11 (Majumdar and Dhar [20]) There establishes a bijection between Al o,nd 71. 

Proof. First we order all edges incident on each site x G Gf' in some order of preference. For each 
conhguration h G AIl, we consider a following discrete-time growth process of graph on Gl, which 
is called a burning process on {GL,h). Let Vq = Vq = {r}, Eq = it) and Tq = {Vo,Eq). Assume that 
we have nonempty sets Tt = {Vt,Et) and Vt with t G Nq. Let 


Vt+i = 

For each y G Vt+i, consider 


yeG^l\Vf.h{y)> Y. (-AL(x,y)) 


x:xeGl’lvt 


We must have 


Et+i{y) = |e G G^^ : e connects y and a site in Vt 


h{y)< Y (-Ai,(x,y))|.Et+i(y)|, 
x:xeGl\yi 


}■ 
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since h € Sl. If |£'t+i(y)| = 1, then name that edge as e(y). If |-Et+i(y)| > 2, then write 

Hy)= 

x:xeG^”)\Vt 

and choose the (s + l)-th edge in Et+i{y) as e(y). We define 

Vt+i = VtU I^+i, Et+i = EtU {e(y) : y G and T^+i = (V^+i) -E't+i)- 

By the assumption h G Al, there is a finite time a < oo such that and -Bo- = 

By the construction, To- = (14-,-Bo-) is a spanning tree on Gl- Since this growth process of Tt,t G 
{0,1,-- - ,cj} is deterministic for a given configuration h G Al, it gives an injection from Al to 
7 l. This fact and Lemma 2.10 give 1^1^! < \Tl\ = det A/,. On the other hand, Lemmas 

2.3 and 2.7 give det < \Al\- Then we can conclude \Al\ = det A^ and the 

burning process gives a bijection between Al and 71. | 

Combining Lemmas 2.3, 2.7, 2.10, and 2.11, we have the following proposition. 

Proposition 2.12 For the DASM on Al, TZl = Al- 

3 Avalanche Propagators 

3.1 Integral expressions for propagators 

Define 

GL(x,y) = EL[T(x,y,B)], x,yGA, 

where T(x,y,/i) is given by (1.2) and the expectation is taken over configurations {h} in the 
stationary distribution Pl- Gl(x, y) is regarded as the avalanche propagator from x to y [6]. 
Sometime in an avalanche caused by a deposit of a grain of sand at x, this site x topples many 
times. The set of topplings between the first and the second toppling at x is called the first wave 
of toppling. There can occur many waves in one avalanche and Gl(x,x) gives the average number 
of waves of topplings in an avalanche [15]. 

Consider the stationary distribution Pl of the DASM. For addition of a particle at any site 
X G Al, the averaged influx of grains of sand into a site z G A^, is given by l(z = x) + 
X/y:y^z Gl(x, y)| Ai(y, z)|n, and the averaged outflux of them out of z by Gl(x, z)Ai(z, z)n using 
the avalanche propagators. In P^, equivalence between influx and outflux must hold at any site 
z G A^. This balance equation is written as 

^ Gl(x, y)Az,(y,z) = il(z = x) Vx,z G A^ 
y-ye^L 

and thus the propagator is given using the inverse matrix of A^,. 

Lemma 3.1 (Dhar [6]) 

GL(x,y) = ;i[A^^](x,y), x,yGAi. (3.1) 
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The matrix can be diagonalized by the Fourier transformation from x = (xi,--- ,Xd) to 
n = (m, • • • ,nd), 

f/t{n,x) = n) = exp ■ n) , 

where x • n = Yl'i=i a-s 


E E Ul (n, x) Al (x, y)U^^ (y, m) 

xixeAi y.yehL 


f 1 ^ 

= 2(i < (1 + a) — - ^ cos I 
I i=l 

= AL(n, m), n,mGAL. 


2L + 1 


rii ) > l(n = m) 


Then, (3.1) is obtained as 


GL(x,y) = - ^ 17^^(x,n)[A^^](n,m)t/L(m,y) 


n:n£AL mrmEAi, 


1 1 g-27rV^(x-y)-n/(2L+l) 

= —ww-TT > - 3 -;;- .nonumber (3.2) 

2dn{2L + ir^^^^Jl + a)-{l/d)j:U<,o.{.^n.) 

Lemma 3.2 There exists a limit G(x — y) = lim£,^oo Gl{^, y), x, y G and 




3 — \/^X -0 


2dn f \ J_^ 2 tt (1 + a) - (1/d) cos 9., 
Proof. Consider the Euler-Maclaurin formula for / G C^(M), 


■, xGZ*^. 


^ 1 rb+Mc 1 ..Ml 

Y fib + nc) = - / f{e)de + -[fib) + f{b + Me)] + —Y f^^^ib + c(f^ + </>)), (3-4) 

n=0 '^46 ^ ^^0 

where M G N, 6, c G M, /'■^^(d) is the second derivative of /(0), and 0 < (/> < 1 (see, for instance. 
Appendix D in [1]). Assume that 


fiO) = 




(1 + a) — (l/(i)(cos 9 + 02 ) ’ 

where a, ai, 02 are constants. Applying the Euler-Maclaurin formula (3.4) with b = —27 tL/(2L+1), 
M = 2L and c = 27rj{2L + 1), we have 

^ -27rC=TQ:i(n-L)/(2Li + l) 


^ (1 + a) - (l/d){cos(27;^(n - L)) a2} 


= (2L + 1) 


r2nL/{2L+l) 

I- 2 ttL/{ 2 L+i) 2vr (1 + a) — (1/d) (cos 6* -|- 02 ) 


/ - 


2L + 1 


2L + 1 


1 / 27r 


12 V 2 L-M 


2 2L-1 
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By dividing the both sides of the equality by 2L + 1 and take the limit L '[ oo, we obtain 


lim 




^—2n\/—lain/ (2Li+l) 


Ltoo 2 L + 1 (1 + a) - (l/(i){cos( 2 ^n) + 02 } 


(le 




2 tt (l + a) - (l/(i)(cos0 + a 2 )' 

Repeating this procedure d times, we can prove Lemma 3.2. | 


3.2 Long-distance asymptotics 

Now we consider the asymptotic form in |x| f 00 of G'(x). Here we follow the calculation found 
in Section XII.4 of [21] for the asymptotic expansion of two-point spin correlation function of the 
two-dimensional Ising model. By using the identity 



1 

a 


and the definition of the modified Bessel function of the first kind 


In{z) 




\/—ln(p-\-z cos (f) 


we have 

The asymptotic expansion of In{z) for large n is found on p.86 in [9], 


, , 1 exp [(n^-|-— nsinh ^inlz) , , , 

Uz) = X (1 + 0{l/n)) , 




{v? -i- 


and we obtain 


G(x) = 




d/2 




2dn\27rJ Jq + {s/ dyy/* 

X -|-C)(max{l/xj}) 


exp[- 5 (x,s)] 


(3.5) 


where 


5 (x,s) = (1 a)s - ^ 


i=l 




- ,,-l (i 

+ y Xi smh I -x, 
i=i 


We can evaluate (3.5) by the saddle-point method and obtain the following result. 
Theorem 3.3 Let 

ci(d, a) = 


1 


47r(a -|- 1) 


a{a + 2)d 
27r(a -|- 1) 


1 {rf-3)/2 


(3.6) 
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and 


1 


i{d,a) = 


\/dsinh ^ a(a + 2) 
r/ien, /or i/ie DASM with d > 2, m, n G N, a = ml{2dn), 

j2 


lim-log 

r^oo T 


ci{d, a) 


-G(x(r)) 


/(d,a)’ 


where 


X r = 


= r- 


r>o. 


(3.7) 


(3.8) 


(3.9) 


VVd’ ’Vdy 

Proof. Let 5 ) and 5 ) be the first and second derivatives of s) with respect to s 


1 ^ 


2=1 


1 + I —Xi 
s 


1 1/2 


d 


5^2) (x,s) = 3^x2 


2=1 


d 


1 + -Xi 

s 


- 1/2 


For each x, let so(x) be the saddle point at which g^^\x,s) vanishes, 

5((^)(x,so(x)) = 0. 


(3.10) 


Then 


G(x) = 


1 /1 


2dn \2 ttJ [xf + sq{x)‘^/ d?Y/^ 


exp[-c/(x,so(x))] 


X / du exp 


^ n 


1 


--g^‘^Yx,so{x))u^\ X ( 1 + C)(max{l/xi}) 


1 


2dn \27rJ (x? + so(ai)^/<^^)^'^"^ 

271 ^ 

ff(2)(x,So(x)) 


exp[- 5 (x,so(x))] 


X ' 1 + 0(max{l/xi}) . 


Here we can prove that the higher derivatives of g{x, s) only give the contributions of order 
0(maxj{l/xi}). See p.304 in [21]. Now we consider the case 

in which e/s are finite and fixed and r 3> 1. The equation (3.10) for the saddle point is now 






2=1 


So(x)2 


. 1/2 

A £i] I = (1 + a)d, 
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and it is solved as 




This gives 


9 (x.,„W) = + + 


Vdr smh ^ \/ a{a + 2) + sinh ^ \/a(a + 2) x + 0{l/r) 


<,« (X, ,„(x)) = 1 +o(i/,2) 

y/d a + 1 r 


Then we have the estimation 


ci(d,a) 1 r 


d 

—- - X{a)y^ Ei X (1 + C)(l/r)), as r t oo 

^ J 


for X = {r/y/d + El, ■ ■ ■ ,r/Vd + Ed), where ci{d,a) and ^{d,a) are given by (3.6) and (3.7), respec¬ 
tively, and 


A(a) = 


^{d,a) 

= sinh“^ y/a{a + 2) 

= log(l + a + \/ a{a + 2)). 


(3.11) 


If we put Ei = Q,1 < i < d, then G(x) is reduced to be 


G{x{r)) = G{r) x (1 -|- 0{l/r )), as r t oo 


G(r) = 


ci(d,a) e-GgKG 
n rG~^)G 


(3.12) 


It proves the theorem. 


4 Height-0 Density and Height-(0,0) Correlations 


^ r 1 2 , 1 

a,l3 ^ <Q, ,hc - 

n n n 


define 


-Pq:,l(x) = El[1(/i(x) = a)], 

-Pa/3, L(x,y) = Ei[l(h(x) = a)l(h(y) =/?)], x,y e Al- 
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Pa,L{^) is the probability that the site x has the height an measured in the unit of grain of sand, 
1 /n, and Pa/3,L{^,y) is the {a, j3)-height correlation function [19, 5, 23]. 

For the two-dimensional BTW model on with open boundary condition, Majumdar and 
Dhar [19] proved the existence of the infinite-volume limits 

Po = limPo,L(x), 

L^oo 

-Poo(x(r)) = lim Poo,L(0,x(r)), 


where x(r) = [r j\f2,r j They gave an 8 x 8 matrix Mi{r), whose elements depend on L and 
r, such that 


^oo,L(0,x(r)) = det ML(r), VL > 


v/2’ 


and showed that every elements converge in the infinite-volume limit L f oo with a finite r. Then 
the matrix M{r) = ^L{r) is well-defined and we have the determinantal expression 

Poo(x(r)) = det M(r). 


Moreover, they showed that 

limPoo(x(r)) = 
rtoo 

and 

Coo(x(r)) = ~ cx -Ir-", as r t oo. (4.2) 

Majumdar and Dhar claimed [19] that the result (4.2) is generalized for the d-dimensional BTW 
model with d > 2 as 

C'oo(x(r)) ~ as r t oo. (4.3) 

In an earlier paper [28], all these facts also hold for the two-dimensional DASM, if we prepare 
10 X 10 matrix Mi{r). (See also [5] and [23] for other generalizations of [19].) Here we show the 
result for the height-0 density and the height-(0,0) correlations of the DASM with general d > 2. 


4.1 Nearest-neighbor correlations 

First we prove the following Lemma. 

Lemma 4.1 Any configuration h G Sl, in which there are two adjacent sites zi, Z 2 G A/,, |zi—Z 2 
1 , such that h{zi) < 1 and h{z 2 ) < 1, is not allowed. 

Proof. Let F = {zi,Z 2 } C Al. Then 

^ (-Ai(x,Zi)) = -Al(z 2 ,Zi) = 1, 

x:x£F,x^zi 

and 

^ (-Ai(x,Z 2 )) =-Al(zi,Z 2 ) = 1, 

x:xeF,X7^Z2 

by (1.1). Then if h{zi) < 1 and h{z 2 ) < 1, the condition of FSC (2.12) is satisfied. | 

By Propositions 2.6 and 2.12, the above lemma implies the following. 


16 



Proposition 4.2 For any L >2, 


-Pa/3,L(0,±ei) = 0, l<i<d, a ,/3 G jo, -,-,...,1- -1. 

I n n n) 


Then, 


Pai3{0, ±ej) = lim Pa/3, l(0, ±ej) = 0, 1 < i < d, a, /? G jo, 

L^oo 


1 2 1 

n n n 


4.2 Determinatal expressions of Po,l(0) and Foo,l(0,x) 

Let Oj, 1 < i < d be the i-th unit vector in Define a real symmetric matrix with size (2L + 1)*^ 


as 


pf(v,w) = <: 


—Zic + l/n, if V = w = 0, 

-1, if V =’w, |v| = l,v /-erf, 

-1 + 1/n, if V = w =-erf, 

1, if V = 0, |w| = 1, w 7 ^-erf, 

1 — 1/n, if v = 0,w = —erf, 

0, otherwise, 

where v, w G A/,. 

Lemma 4.3 Let El he the unit matrix with size {2L + l)*^. Then 

Po,l( 0) = det (Pl + nGiPf) . 

Proof. Define a set of allowed configurations conditioned h{0) =0, 

= {h£AL: h{0) = 0}. 

By definition (4.1), Proposition 2.6 with Lemma 2.3 and Proposition 2.12 gives 


Po,l(0) = 


iA“’i 




n 


(4.4) 


Assume that h G . Then as shown in the proof of Lemma 2.11 we can uniquely define a burning 
process Tt, t G {0,1,... ,^ a} on {Gl, h) associated that Tj becomes a spanning tree on Gl at time 
t = a. Define a configuration h' as 

{ /i(z) - 1, if |z| = l,z/-erf, 

h{z) - 1 + 1/n, if z = -erf, 

/i(z), otherwise 

for z G Al- Now we consider a new DASM which is defined by the matrix Af given by 

A'l = Al + B^l\ (4.5) 
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and let be a set of all allowed configurations of this DASM and G'^ be an associated graph 
to (Ai, A'^). Then we consider a burning process T/ = (F/, E^),t G {0,1,... , ct} on h'). By 
dehnition of and h!, we can make 

Vt = Vl VtG{0,l,...,a}, 


and gives a spanning tree on G'j^. By Lemma 2.11, this means h' G A'j^. Since there is a bijection 
between h and its associated burning process Tt, t G {0,1,..., ct}, we have a bijection between A® 
and A']^. By Lemmas 2.10 and 2.11, |A®| = \A!i\ = det A'^. Combining (4.4) and (4.5) 

gives 


-Po,l(0) 


det 

det A^, 
det(A^iAl) 
det(L;L + A-iB®). 


Then we use Lemma 3.1 and the proof is completed. | 

Next we consider the two-point function i-oo,L(0,x), where we assume that 2 < |x| < L. We 
define a real symmetric matrix with size (2L + 1)*^ as follows. For v,w G A^, 


B 


(O.x) 


(v,w) = 


-he + 1/n, 

if 

v = w = 0orifv = w = x. 

-1, 

if 

V = w, |v| = 1, V 7 ^ -ea, 


or 

if V = w, V — X = 1, V 7 ^ X — e^, 

-1 + 1/n, 

if 

V = w = —e^, or if v = w = x 

1 , 

if 

V = 0, |w| = 1,’w 7 ^ -erf. 


or 

if V = X, w — X = 1, w 7 ^ X — Orf 

1 - 1/n, 

if 

V = 0,w = -erf. 


or 

if V = X, w = X — erf. 

0, 

otherwise. 


Following the same argument as Fb,L(0) we can prove the next lemma. (See Fig.6.) 


Lemma 4.4 For 2 < |x| < L, 

Poo,l{0,^) = det [El + nGiSf. 


4.3 Infinite-volume limit 

Since the number of nonzero elements of (resp. ) is only Gd -|- 1 (resp. 2(6d -|- 1)), we 

can replace the matrix El + riGiB^^^ (resp. El + nGiB^^’^^) with size (2L + l)'^ by a matrix with 
size {2d + 1) (resp. 2(2d -|- 1)) without changing the value of determinant. Explicit expressions are 
given as follows. 

Let 

to, if f = 1, 

q* = < ej_i, if 2 < i < d-k 1, 

[ if d + 2<i<2d-kl. 
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Figure 6: The matrix A'[ = is considered for Poo,l(0,x) with |x| = r. In the corre¬ 

sponding graph G'l the site 0 (resp. x) is connected to —(resp. x — e^) by a single edge, but all 
other edges between 0 (resp. x) and its nearest-neighbor sites are deleted. 


Define a matrix )i<i,j< 2 d+i with elements 

= GL(0,x-Fqj -qj), l<i,j<2d+l. (4.6) 

We also define a real symmetric matrix B = (i3ij)i<ij<2d+i with elements 




-hc + l/n, 
- 1 , 

-1-bl/n, 

1 , 

1 - 1/n, 

0 , 


if i = j = 1, 
if 2 < i = j < 2d, 
if i = j = 2d + 1, 
i = l,2<j <2d, 
if i = l,j = 2d + 1, 
otherwise. 


Then define 2(2d -I- 1) x 2{2d + 1) matrices 


e(^^(0,x) 


^(D(0) g{D(x) \ 

ig(-^)(x) g(^)(0) J ’ 


X E Al, 


where is a transpose of ^(^^(x), and 


B 


P 0 \ 

0 B )' 


We have 

Po,l(0) = det [e + (4.7) 

and 

^oo,l(0,x) =det (p + ng(^)(0,x)p) , (4.8) 

where E denotes the unit matrix with size 2d+l in (4.7) and with size 2{2d+l) in (4.8), respectively. 

It should be remarked that the sizes of the matrices in the RHS’s are independent of the lattice 
size L and determined only by the dimension d of lattice. The dependence of L is introduced 
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only through each elements of Q^^\x.) given by (4.6). Lemma 3.2 guarantees the existence of 
infinite-volume limit L f oo of these elements and we put 

Gij(x) = lim Qi^\x) = G(x Qj- - qi), I < i,j < 2d + 1, 

G{x) = (0ijr(x))i<jj<2rf+l, 

6(0,X) = li.n6™(0,x)=(fj")j 

where G(x) is explicitly given by (3.3). Then we have the following. 

Proposition 4.5 There exist the infinite-volume limits 


Po = lim Po,l(0), Poo(x) = lim Poo,l(0,x), x G IT, 
Z,t°o Lt°o 


and they are given by 
and 


Pq = det {E + nQ{0)B) 
Poo(x) = det (^E + nG{0, x)BJ , 


X G Z"*. 


4.4 Evaluations of determinantal expressions 

From the determinantal expressions of Pq and Poo(x) given in Proposition 4.5, the following explicit 
evaluations of these quantities are obtained. 


Theorem 4.6 (i) Define 


71 = 


d 


dOi 


2d 


J-^ 27r (1 + a) - ( 1 /d) Ecos ( 


and 




dOi 


,-2 v^(01 +02) 


2d l\ 27r (1 + a) - ( 1 /d) Zti cos 9,' 
Then, for the DASM with d > 2, m, n G N, 

1 — 2da7i 


Pn = 


2 dn 


[2{1 - d( 7 i - 72 )} + (1 - 4d7i)a - 2 d 7 ia^ 


X [2(d - l)( 7 i - 72 ) - (1 - 4d7i)a -|- 2d7ia^] 

X [{l-(7l-72)}'-{(2d(l + a)2-l)7l-(2d-l)72-(l + a)}2]"-^ (4.9) 

where a = m/{2dn). 

(ii) Let 


C„„(x) = 


p2 


X G Z“'. 


(4.10) 
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Then, there exists a nonzero factor C 2 {d, a, n) such that for the DASM with d > 2, m, n G N 


lim — log 

r^oo T 


r^d 1 

C2{d,a,n) 


C'oo(x(r)) 


2 

fid, a)' 


(4.11) 


where a = m/{2dn), fid, a) and x(r) are given by (3.7) and (3.9), respectively, and that 


where 


C2id,a,m/i2da)) _ / d d{l + id - 1)^} ^ 

04,0 a(<^+i)/2 y 27 r 2 y 27 r((i — 1)7 



d9i 1 - e“2v^('?i+®2) 
27r 1 - (l/d) Y.i=i cos^i 


(4.12) 


In the following, we will explain how to prove this theorem. Let 

= E + nQiO,xir))B, r > 0, x(r) G 


where is a unit matrix with size 2(2(i + 1). That is, 



rh^^\r) \ 

mA) J ’ 


where for 1 < i < 2d + 1 



l(i = 1 ) + EHV "fe( 0 ) 

-{(1 - 1/n) + hcjnQiiiO) - 0 i 2 d+iiO), 
l(i = j) + nlQiiiO) - Qij iO)], 

, l(z = 2d + 1) + (1 - l/n)nigniO) - gi2d+ii0)j, 


if J = 1, 
if 2<j< 2d, 
if j = 2d + 1, 



= n X < 


EiiV&wo) 

-{(1 - 1/n) + /ic}0ii(x(r)) - (l/n) 0 i 2 d+i(x(r)), 
^ii(x(r)) - aij(x(r)), 

(1 - l/n)(aii(x(r)) - t/i2d+i(x(r))). 


if j = 1, 
if 2 < j < 2(i, 
if j = 2d + 1, 



= n X < 


ES'5H(x(r)) 

-{(1 - 1/n) + 7 c} 0 ii(x(r)) - (l/n)e 2 d+ii(x(r)), 
0ii(x(r)) - aji(x(r)), 

(1 - l/n)(aii(x(r)) - 02^+1 i(x(r))). 


if / = 1, 
if 2 < j < 2d, 
if j = 2d + 1. 


We find that 

2d+l 

^ = 1 - 2dan0ii(O), 

1=2 

2d+l 

^!i^(^) + X] = -2dant/ii(x(r)), 

1=2 

2d+l 

^iX^) + X/ ~ —2dan^ii(x(r)), 1 < i < 2d + 1. 
1=2 
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For 1 < i < 2(i + 1, let 


niij = 


1 — 2danQii{Q), if j = 1, 
( 1 ) 


m. 




^ j -2danGii{yi{r)), if j = 1, 

= 1 <■>(.), 


I -2(ian0ii(x(r)), if j = 1, 

= 1 4'F), 


if 2 < J < 2d + 1, 

if J = 1, 
if 2<j<2d + l, 

if j = 1, 
if 2<j<2d + l. 


Then 


Pq = det = det m, 

Poo(x(r)) = det (r) = det M(r) with M{r) = 

Note that, if we introduce the the dipole potential 

hh,ji),{i2,j2)i^i'^)) = ^hii(x(r)) -Gi2j2i^ir)), 1 < ii,i2,ji,j2 <2d+l, 

the elements of the matrix M{r) are expressed as follows; for 1 < i < 2d + 1, 


nii 


1 — 2dan0ii(O), 

= { i(^ = i)+ ^</'(i,i),(i,i)(o)> 


if j = 1, 
if 2 < j < 2d, 

l(i = 2d+1) + (1 - l/n)n(/)(i^i)^(i_2d+i)(0), if j = 2d+l, 


(4.13) 


(4.14) 


mij{r) = n X 


-2da0ji(x(r)), 

(1 - l/n)(()(i_i)^(j^2d+i)(x(r)), 


if j = 1, 
if 2 < j < 2d, 
if j = 2d + 1, 


rhij{r) = n X 


-2da0ij(x(r)), 

<('(i,i),(i,i)(x(r')), 

(1 - l/n)(/>(i^j)_(2d+i,i)(x(r)). 


if j = 1, 
if 2<j< 2d, 
if j = 2d + 1. 


Now we study the asymptotics of Fbo(^) in r f oo. Theorem 3.3 and its proof given in Section 
3 implies that with any finite Cj’s, 


G 



G(r)exp —A(a)^^( 


2 = 1 


x(l + 0(l/r)). 


as r t oo 


with (3.6),(3.7), (3.11), and (3.12). Then we see 


rh(r) = nG{r)n{r,\){l + 0{l/r)), 

m{r) = nG{r)n{r,—X){1 + 0{l/r)), as r t oo. 
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where n{r,X) = (njj(r, A))i<ij< 2 d+i with elements, 


nij{r,X) 


—2da, 

(1-e-^), 

(1 -e^), 

(l-l/n)(l-e^), 

—2dae^, 

—2dae~^, 
e'^(l — e""^), 
e^(l — e'*'), 

(1 — l/n)e^(l — e^), 
e“'’'(l — e"'*'), 
e“'’'(l — e"^), 

(1 — l/n)e“'’'(l — e'*'), 


if i = j = l, 

if i = 1,2 < j < d + 1, 

if i = l,d + 2 < j < 2d, 

if i = l,j = 2d + 1, 

if 2<i<d+l,j = l, 

if d + 2<i<2d+l,j = l, 

if 2 < i,j < d + 1, 

if 2<i<d+l,d + 2<j< 2d, 

if 2<i<d+l,j = 2d + l, 

if d + 2 <i <2d+l,2 < j < d+1, 

if d + 2 < i,j < 2d, 

if d + 2 < i < 2d + l,j = 2d + 1. 


We obtain a matrix M'{r) from M(r) by subtracting (the first row) xe"^ from the i-th row with 
2 < i < d+1, (the first row) xe"'*' from the i-th row with d + 2 < i <2d+l, (the {2d + 2)-th row) 
xe”"^ from the i-th row with 2d + 3 < i < 3d + 2, and (the {2d + 2)-th row) xe'*' from the i-th row 
with 3d + 3 < i < 2{2d + 1). We have 


with 


m’ij{X) 


M'{r) 


m'{X) fh!{r,X) \ 
m!{r,—X) m'{—X) J 


1 — 2danGii{0), 

(1 - l/n)n(/)(i^i)_(i^2rf+i)(0), 

(1 - e^) - 2dan{gii{0) - e^0ii(O)), 

(1 - e“^) - 2dan{gii{0) - e'-^^nl^O)), 

< l(i = j) + n[0(i^i)^(ij)(O) - e^0(i_i)^(ij)(O)], 
l(i = j) + n[0(i^i)^(ij)(O) - e-'^())(i^i)^(ij)(0)], 

(1 - 1/n) 

^^['/’(i,l),(i,2d+l)(0) — e'’'(/)(i^l)^(1^2d+l)(0)]) 

l{i = 2d+ !) + {!- 1/n) 

X^[</’(i,l),(j,2d+l)(0) — e“^(/)(l^l)^(1^2d+l)(0)]) 


if 

if 

if 

if 

if 

if 

if 

if 

if 


i= j = 1, 
i = l,2<j <2d, 
i = l,j = 2(i+ 1, 

2<i<d+l,j = 1, 
d + 2<i<2d + l,j = 1, 
2<i<d+l,2<j <2d, 
d + 2 <i <2d + l,2 < j <2d, 

2<i<d+l,j = 2d + l, 

d + 2 <i <2d + l,j = 2d + l, 

(4.15) 


and with 


fh[j{r. A) = nG{r) x < 


-2da{l + 0{l/r)), 

(l_e-A)(i + 0(l/r)), 
(l-e^)(l + 0(l/r)), 

(1 - l/n)(l -e^)(l +0(l/r)) 
0{l/r), 


if i = j = l, 
if i = 1,2 < j < d + 1, 
if i = l,d + 2 < j < 2d, 
if i = l,j = 2d + 1, 
otherwise, 


(4.16) 


so that 

PoO(x(r)) = det M(r) = det M'(r), r > 0, x(r) G 
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Now we expand det M'{r) along the first and the (2d+2)-th rows. Let A:)| be the determinant 

of M'{r) with the first and the {2d + 2)-th rows and the j-th and the k-ih. columns removed and 
multiplied by — (—x (|_i) 2 d+ 2 +A: _ Then we have 

2(2d+l) 2(2d+l) 

det M'(r) = E E M' {r)ijM' {r)2d+2,k\M' {j,k)\. 

j=l k=l,k^j 


Remark that, by (4.15) and (4.16), 

\M'{j,k)\ = 0{l/r), as r —)• oo, 


if 1 < j. A: < 2(i + 1 or 2(i + 2 < /c < 2{2d + 1), and 

\M'{j,k)\ = |m'(^')(A)| X |m'W(A)| x (l + 0(l/r)), as r —)■ oo, 

if 1 < i < 2(i + 1 < A < 2{2d + 1) or l<k< 2d+ l<j< 2{2d + 1), where is the 

(1, j)-cofactor of m'{X). Then 


det M'{r) 


2d+l \ / 2d+l 

m'y(A)|m'(^)(A)| Y 
i=i / \ i=i 


( 2d+l \ / 2d+l 

Y^ 'm[jir,X)\m'^^\-X)\\ ( ^ my(r,-A)|m'(^)(-A)| 

= detm'(A) x detm'(—A) + det?fi(A) x det?77.(—A) x (nG'(r))^ (1 + 0{l/r)), (4.17) 


where m{X) = (mjj(A))i<ij< 2 d+i with elements 


mij{X) 


' -2da, 

1 — e^, 

1 — e~'^, 

(1 - l/n)(l -e"^), 

. m'ijiX), 


if i = j = l, 
if i = 1, 2 < j < d + 1, 
if i = l,d + 2<j< 2d, 
if i = 1, j = 2d + 1, 
otherwise. 


We find that 


detm'(A) = detm'(—A) = detm. 


The determinantal expressions (4.13) with (3.12), (4.17), and (4.18) give 


lim Poo(x(r)) 

r^oo 


lim{(det m)^ + det m(A) det m(—A)(nG(r))^} 
(det m)^ = Pq . 


Here we set 


detm(A) = a detm* (A), 


(4.18) 
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with a matrix m*{X) = {rn*j{X))i<ij< 2 d+i with elements 


m*j{X) = { 


-2d, 

(1 — e'’')/a^/^ 
(1 




if i = j = 1, 
if i = 1,2 < j < d + 1, 
if i = 1, d + 2 < j < 2d, 
if i = 1, j = 2d + 1, 


(1 —l/n)(l —e 

(1 — 

-2da^/‘^n{Qii{f)) - e^0ii(O)), 

(1 — e“'^)/a^/^ 

-2da}/'^n{Qii{f)) - e“^0ii(O)), 
l{i = j) + 

l{i = j) +ni(/>(i_i)_(jj)(0) - e-^0(i_i)^(ij)(O)], if d + 2 < i < 2(i + 1, 2 < j < 2(i, 
(1 - 1/n) 

l(i = 2d+l) + (l-l/n), 


if 2 < i < (i + 1, j = 1, 

if(i + 2<i<2(i + l,j = 1, 
if 2 < i < d + 1,2 < j < 2(i, 


(4.19) 


if 2 < i < d + 1, j = 2d + 1, 
ifd + 2<i<2d + l,j = 2d+l. 


By the definition (4.10), we see 


C'oo(x(r)) = a 


2 det m*{X) det m*(—A) 


{nG{r))‘^ X {1 + 0{l/r)), as r too. 


(det m)2 

Since G{r) is given by (3.12), (4.11) of Theorem 4.6 (ii) is proved with 

2 det m* (A) x detm*(—A) 


C 2 {d,a,n) = {aci{d,a)y 


(det m) 


Now the problem is reduced to the calculation of detm and detm*(A). Consider a matrix 
R = {Rij)i<i,j<N with elements 


Rij = < 


u, 

b, 

c, 

(1 - l/n)c, 

q, 

e, 

f, 

l + v, 
h, 

(1 - 1/n)/, 

(1 - l/n)h, 

s, 

t, 

1 k, 

(1 - l/n)s, 

1 + (1 - l/n)k, 


if i = j = 1, 

if i = 1, 2 < J < d + 1, 

if i = 1, d + 2 < J < 2d, 

if i = 1,/ = 2d + 1, 

if 2 < i < d + 1, j = 1, 

if d + 2 < i < 2d + 1, j = 1, 

if 2<i<d+l,2</< 2d, j ^ i,j ^ i + d, 

if 2 < i = / < d + 1, 

if 2<i<d,j = i + d, 

if 2 < i < d, j = 2d + 1, 

if i = d + 1, / = 2d + 1, 

if d + 2<i<2d+l,2</< 2d, j y i,j y i — d, 

if d + 2 < i < 2d + 1, j = i — d, 

if d + 2 < i = / < 2d, 

if d + 2 < i < 2d, j = 2d + 1, 

if i = j = 2d + 1. 


(4.20) 


We perform the following procedure on R. 
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(i) Subtract (the first row) xq/u from the f-th row with 2 < i < d + 1. 

(ii) Subtract (the first row) xe/u from the i-th row with d + 2<i<2d + l. 

(iii) Subtract the second row from the i-th row with 3 < f < d + 1. 

(iv) Subtract the (d + 2)-th row from the f-th row with d-|-3<f<2d-|-l. 

(v) Add the j-th column to the second column with 3 < j < d + 1. 

(vi) Add the j-th column to the (d -|- 2)-th column with d -|- 3 < j < 2d. 

(vii) Add (the (2d -|- l)-th column) xl/(l — 1/n) to the (d -|- 2)-th column. 

(viii) Subtract (the (d -|-j)-th column) x(f — s)/(l -|- A: — s) from the j-th column with 3 < j < d. 

After these procedures, by changing the orders of rows and columns appropriately, we obtain the 
following identity. 

t — 1 

det R = u X 1 + V — f — - - - - {h — f) x (1 -|- /c — sY~‘^ x det S, (4-21) 

L K — S 

where S = (5'ij)i<ij<4 with elements 


5ii = 1 -t- u -|- (d — 1)/ — dbq/u, 

-5i 2 = h -1- (d - 1)/ - dcq/u, 

5i 3 = (1 - l/n)(/ - cq/u), 

-5i4 = / - bq/u, 

S 21 = t + {d — l)s — dbe/u, 

S 22 = 1 + k + {d — l)s — dceju. 

S 23 = (1 - l/n)(s - ce/u), 

1 

II 

S 31 = 0, 

S 32 = l/{n - 1), 

-533 = 1 + (1 - l/n){k - s), 

1 

II 

S 41 = 0, 

to 

II 

0 

-543 = (1 - l/n)(/i - /), 

844 = 1 + v - f. 


Define 

gQ = nG{Q), gi = nG{ei), g2 = nG{2ei), g^ = nG{ei + 62), 

where G'(x) is given by (3.3) and ei, 62 are the unit vectors in the first and second directions in 
Since the system is isotropic, we can find that the matrix m defined by (4.14) is in the form (4.20) 
with 

u = 1 — 2daqn, h = c = qn — qi, 

q = e = l- 2dagi, f = s = gi - gs, (4.22) 

v = k = gi-go, h = t = gi-g 2. 

By Lemma 3.2 and the isotropy of the system gives 

2d(l a)go - 2dgi = 1, 

2d(l 0)5-1 - (5-0 + g2 + 2(d - 1)5-3) = 0, 

which are written as 

/ N 1 

5-1 = {l + a)go-^, 

52 = [2d(l + a)2-l]5o-2(d-l)53-(l + a). (4.23) 
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The formula (4.21) with (4.22) and (4.23) gives 

Po = detm = [2{l - d{gQ - g^i)} + [I - MgQ)a - 2dgQa] 

X [2{d - l)(5(o - 53 ) - (1 - 4(i5fo)a + 2dc/oa^]^ 

X [{1 - (50 - S's)}^ - (52 - S's)^]'^ ^ • (4.24) 


It proves (4.9) of Theorem 4.6 (i). 

It should be noted that, if we put n = 1 and take a 0 limit in (4.24), we have the formula 
Po = ^ ^ ^ (1 - d503)503[(l - 503 )^ - 523]'^~^> 

where 

503 = lim(5o - 93 ), 923 = lim(52 - 93 )- 

a 4-0 04,0 

In particular, 503 = 1/vr and 523 = 1 — l/'7rfor(i = 2 [27], and thus we have 



This coincides with the value of Pq obtained by Majumdar and Dhar [19] for the two-dimensional 
BTW model. 

We can also find that the matrix m*{X) defined by (4.19) is in the form (4.20) with 


u = —2d, 

6 = (1 — e^)/a^/^. 

c = (1 — e~^)/a}P, 

5 = (1 — e^)/a}P — 2da}P{g 

e = (1 - e~^)/a}P - 2da}P{gi - e"^5o), 

f = (di - 93 ) - e^igo - 9i), 

s = {91 -53) -e“^(5o -di), 

II 

1 

0 

1 

CD 

0 

1 

k = {gi- 50) - e"^(5o - di), 

II 

1 

to 

1 

CD 

> 

0 

1 

II 

1 

to 

1 

(TD 

1 

> 

0 

1 


The formula (4.21) gives 



detm*(A) = -2d [{1 - (50 - 5 - 3 )}^ - {92 - 53 )^]'^ ^ x det S’, 


where ^ 

det S = bi(d, a, A) + b 2 (d, a. A) — . 

n 

with some functions bi and 62 of d,a,X. Since (3.11) gives 

gA(a) _ Q _j_ y/ a{a + 2) = 1 + V2a^P + 0 {a), as a I 0, 

we found that 


6i(d,a,A) = 0(0^), 

b2id,a,X) = ^(50 - 5'3){1 - 450-5'3)}{1 + (d-l)(5'o - 53 )}+ C>(a^^^), as a| 0. 

d 

Thus we obtain 

det m*(A) det m*(—A) r2(i{l-|-(d — 1)503}] ^ 

hm-—--^5-= -—-r-^- 

a4,o (detm)^ [ (« —1)503 

Since lim^^o ci(d, = (d/(27r^))(‘^“^)/"^/(47r), (4.12) of Theorem 4.6 is proved. 
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5 Discussions 


5.1 Critical exponent Ua 

The results (3.8) of Theorem 3.3 and (4.11) of Theorem 4.6 mean that both of G(x(r)) and C'oo(x(r)) 
decay exponentially as increasing r with a correlation length ^{d,a). Since ^(d, a) < oo for any 
a > 0, the stationary state of the DASM is non-critical [28]. Moreover the theorems imply that, if 
we make the parameter n be large with a fixed m, then the value of a = m/{2dn) can be small and 

p—r/^{d,a) 

nG(x(r)) ~ ci(d)a^'^~^^/^ ^(d-i )/2 ’ 

„- 2 r/i{d,a) 

C'oo(x(r)) ~ C 2 (d)a(‘^+^)/^——, as r f oo, (5.2) 

where ci(d) = (d/(27r^))(‘^“^)/^/(47r) and C 2 (d) is given by (4.12). 

Consider a series of DASMs with increasing n with a hxed m. Then we will have an increasing 
series of correlation lengths {^(d, a)} and we will see the asymptotic divergence, 

i{d,a) ~ as a —)• 0 (5-3) 


with 


Va 


— for all d > 2. 
2 


We notice that, if we identify a with a reduced temperature 


t = 


T-T, 

Tc 


(5.4) 


(5.5) 


around a critical temperature Tc in the equilibrium spin system, (5.1) with (5.3) and (5.4) is 
exactly in the Ornstein-Zernike form of correlations in the mean-field theory of equilibrium phase 
transitions (see, for instance, Eq.(61) in Section 3.1 of [14]). This implies that we can regard (5.3) 
as a critical phenomenon with a parameter a approaching to its critical value Uc = 0 and we can say 
that the associated critical exponent Va is exactly determined as (5.4). Vanderzande and Daerden 
discussed the exponent Va for the DASM on more general lattices [29]. 

This exponent may be identified with the critical exponent = 1/2 obtained by Vespignani 
and Zapperi by the generalized mean-field theory [30]. They claimed that they made only use of 
conservation laws to evaluate = 1/2 and thus at least on this result their mean-field theory is 
exact for any d > 2. The present work justifies their conjecture. We can conclude that with respect 
to the avalanche propagators and height-(0,0) correlation functions the upper critical dimension of 
the ASM is two. This result does not contradict to the result by Priezzhev [24], since he studied 
the intersection phenomena of avalanches and for them the upper critical dimension is four. 

The results (5.1) and (5.2) suggest that there exists a scaling limit such that 

lim r'^“^nG(x(r)) = Fg{k), 

r'\oo^a\Q: 

lim r^‘^C'oo(x(r)) = Fc{n), 0 < k < oo 

ri'oOjOj.O: 

(PI'^r=Kly/2d 
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with 

Tg{^) = 

Tc{k) = K'^+^e-^, 

This observation is consistent with the statement 

G(x(r)) ~ as r t oo (5.6) 

and (4.3) claimed by Majumdar and Dhar [19] for the self-organized criticality realized in the d- 
dimensional BTW model with d > 2. (Note that for the two-dimensional BTW model, G(x(r)) — 
G(0) ~ —(l/27r) log r, as r t oo.) 

5.2 The g —)■ 0 limit of the Potts model 

Majumdar and Dhar [20] discussed the relationship between the ASM and the ^ j, 0 limit of the 
q-state Potts model. For g G {2,3, ...}, the g-state Potts model on the lattice Gl = (G^\g^^) 
given by Dehnition 2.8 is dehned as follows. At each vertex v G = AiUjr}, put a spin variable 
s(x) G {1,2,.. . ,g}. The Hamiltonian for the configuration s = {s(v)}^g^(«) is given by 

'H{s) = - l(s(v) = s(w)). 

e={v,w}eG((' 

The partition function of the Potts model in the Gibbs ensemble with a temperature T > 0 is 
dehned by 

Z{q,T) = 

r.i'’) 

s&{l,2,...,q}^L 

= n l + xl(s(v) = s(w)) (5.7) 

s&{l,2...,qfL^ e={v,w}eG){) 

with X = — 1. We consider a subset of denoted by FI C G^\ Each connected component 

in E is called a cluster. Let c{E) be the number of disconnected clusters of E; E = Ei, where 

Ei n Ej = (l},i ^ j. If a vertex v G G^^ is not connected by any edge in E, we write v ^ E. By 
performing binomial expansions and taking the summation over spin conhgurations in (5.7), we 
obtain the Fortuin-Kasteleyn representation of partition function, 

Z{q,T)= Y, (5.8) 

EcG^^'> 

where jEj denotes the number of edges in E. Note that we can regard (5.8) as a function of g G M 
and T > 0. We consider the asymptotics of (5.8) in the limit g {, 0. The dominant terms in this 
limit should be with E such that c{E) = 1 and {v G G^^ : v ^ E'\ = 0 E contains all vertices 
in E is a spanning subgraph of Gl- If we further take the high-temperature limit T f oo 
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X i 0) ''^6 have only spanning subgraphs with a minimal number of edges, which are just the 
spanning trees. Then we have 

lim q~^Z{q,T) = |7l|, 

Tfco 54.0 

where 71 is the collection of all spanning trees on Gl- As shown in Section 2.4, there establishes a 
bijection between 71 and Al (Lemma 2.11) and Al = (Proposition 2.12). (The relation between 
the g' 1 0 limit of the g-state Potts model with finite temperatures and the ASM is discussed in 
Section 7.2 in [7].) The two-dimensional g-state Potts model shows a continuous phase transition 
associated with critical phenomena at a finite temperature 0 < Tc < oo without external magnetic 
field B = 0, when q = 2,3 and 4 [31]. 

Usual critical phenomena of spin models are specified by the behavior of two-point correlation 
functions for the energy density Ge(r, t, b, L) and for the order-parameter density Ga{r, t, b, L). Here 
r denotes the distance of two points, t the reduced temperature (5.5), b the reduced external field 

-^C 

and L the size of the lattice on which the model is defined. It is conjectured in the scaling theory 
that, if L is sufficiently large and we observe the system in the very vicinity of the critical point; 
t <C 1, 6 <C 1, the correlation functions behave as 


G,{r,t,b,L) = L^-^F,(^^,tLy\bLy^) , 

G^{r,t,b,L)=L^^-B^ (^^,tLy\bLy^') , 


(5.9) 


with the scaling exponents x^,Xa, VeiVa-, and the scaling functions Te,T^- If the system is of d- 
dimensional, the hyperscaling relations Xe + yt = d,Xa + Vb = d hold (see, for instance, [13, 14]). 
From the scaling forms (5.9), we expect the power-law behavior of correlation functions at the 
critical point (t = 6 = 0, L t oo) such that 


G,(r) 


— 2 X(: 


G^t) 




as r t CO, 


and in the off-critical regions with L '[ oo, the correlation length ^ = ^{t,b) behaves as 

^(t,0) ~ t~''^ with ut = —, 

yt 

^(0,6)~6“^*’ with = ast|0 ,6|0. 

yb 

For the two-dimensional (/-state Potts model, the critical exponents are determined as functions 
of q through the parameter 

u = u{q) = — cos 

TT 



as [31] 


Xe 

Xa 


1 -|- U 


4(2-n)’ 


yt = ‘2-x^ 
yb = 2-x„ 


3(1 -u) 

2-u ’ 

3(1 — tt)(5 — u) 
4(2 - u) 
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They give the limits 


^ 2, yi —y 0, X/j —y 0, yi) —y 2, cis o' ^ 0 \ y ti ^ 1. 


Majumdar and Dhar [20] noted by their results (4.3) and (5.6) for the BTW models that the 
avalanche propagator G(x(r)) and the height-(0,0) correlation function C'oo(x(r)) in ASM play the 
roles of the order-parameter density correlation function Ga{r) and the energy density correlation 
function Ge{r) in the critical phenomena, respectively. In particular, in the two-dimensional case, 
the power-law exponents are respectively given as 


2Xa 

=0=d-2 

, 2 a;e 

CN 

II 

II 


940 

d=2 

I 94 O 1 


Our interpretation of the present result (5.4) is that introduction of dissipation to the ASM may 
correspond to imposing an external magnetic field B to the Potts models and hence i^a = 1/2 is 
identified with 




940 


1 

Vb 


940 


We remark that the critical exponents for the specific heat a, for the order parameter /3, and 
for the magnetic-field susceptibility 7 of the 


2(1 - 2u) 
~ 3(1 -u) 


—)• — 00 , 


/3 


1 + u 

12 


1 7 — Au + v? , 

R’ = i -asglO^^utl- 

6 6(1 — u) 


We suspect some interpretation of the value = 1/6 in the DASM. 


5.3 Recent topics on height correlations 

In Section 4 the one-point and the two-point correlations of height-0 sites were calculated for 
the DASM with general d > 2. In the two-dimensional case, the three-point and the four-point 
correlations were also calculated for height -0 sites and general property of ‘the height -0 field of 
ASMs’ have been extensively studied from the view point of a c = —2 conformal held theory [18, 8]. 

For the two-dimensional BTW model, in which the values of stable height of sandpile are 
h = 0,1,2, and 3, the height correlations have been calculated also for h > 1. Priezzhev determined 
Pa for a G {0,1,2,3}, where the results with a > 1 are expressed using multivariate integrals of 
determinantal integrands [23]. Poghosyan et al. [22] claimed that the height-0 state is the only one 
showing pure power-law-correlations and that general form of height correlations for h > 1 contains 
logarithmic functions. They showed that for a > 1 

<^ 0 a(x(r)) =-——-~ —(cilogr-Fc 2 ), as r t 00 

with some constants ci,C 2 . Moreover, they predicted that Gay(x(r)) ~ log^r/r^ if a > 1 and 
/3 > 1. These results are discussed with the logarithmic conformal field theory. See also [11]. We 
will see a lot of interesting open problems concerning height correlations for the BTW models and 
the DASMs in higher dimensions. 
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